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ABSTRACT 



The decision as to whether a medical doctor, or other 
expensive specialist, should be carried aboard ship depends upon 
demand for service, consequences of not providing this service, 
and cost of providing the service. We supply a simple preliminary 
mathematical model to aid in making this decision wisely. 
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1 . Introduction . 

The purpose of this paper is to initiate study of the question: 
should expensive specialists, e.g. medical doctors, be assigned to 
ships or submarines on patrol? The issues that arise are as follows. 

(1) If a man is injured or becomes seriously ill while a 
mission is in progress, and if no doctor is present, it may be necessary 
to transfer the individual from the ship to a hospital for treatment. 
Presence of a doctor will, at least in some cases, permit the treatment 
to take place on the ship. Thus the cost of the transfer may be 
eliminated — at, of course, the expense of maintaining the doctor. 

Similarly, but not entirely analogously, breakdown of key 
equipment entails loss of military effectiveness and may require that 
the mission be prematurely terminated. If a skilled specialist is 
aboard such events may be forestalled. 

(2) The cost of retaining doctors or other trained specialists 
is very high. The cost for such a specialist may rightfully include 
some of the expense of his original recruitment and training. 

One subsidiary purpose of this paper is to show that the use 
of simulation methods is not always essential when dealing with problems 
of some complexity, and, in particular, those that involve chance 



effects . 
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2 . The Occurrence of Demands , 

Basic to the question of whether a doctor or highly trained 
and expensive repairman should be added to a ship T s complement is the 
likelihood of demand for services that he alone can supply. We pro- 
pose some probability models for this question. Actually, a complex 
variety of sources may conspire to cause demand. 

Model 1 . Simple Chance Demand, Caused by Accidents or Sudden Disease. 

Imagine that a ship has n individuals aboard when it sets 
out on a mission of duration M. Each individual is thought to have 
a constant probability Adt of experiencing an accident or sudden 
severe illness between t and t + dt, dt being a small number. 

The occurrence of accidents or illness is first assumed to be independ- 
ent from individual to individual. Then our simple model implies that 
each individual experiences his demanding event at time T_^ (i = l,2,...,n) 
measured from the start of the mission, T^ being distributed in 
accordance with an exponential distribution with rate parameter A: 



P{T i <;t} = l- e At 0 £ t < 00 

( 2 . 1 ) 



- 0 



t < 0. 



Next, the occurrence of the smallest T_^ on the ship of crew size 
n is the distribution of the minimum of a sample of n independent 
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T\ ! s, which is exponential with parameter n . Let us call this 

time x ; then 
n 



P{ t > t} 

n 



-nXt 

e 



0 ^ t < 00 



= 1 



t < 0 



( 2 . 2 ) 



According to this model, if no doctor (or repairman) is present: 

(A) The probability that the mission does not terminate during M 
for the cause associated with X is 

P{t > M} = e' nXM (2.3) 

n 

(B) The expected time to the scheduled end of a mission that 
involves a rescue (or terminates early) is 



E[M - t 




T < M] 
n 



rM 

0 



1 - 



(M-x) e 
1 - e 

_M 

-nXM 



-nXx 



-nXM 



_1_ 

nX 



e 



nXdx 



(2.4) 



It is of interest to see what this formula approaches as X becomes 
small, a condition likely to be true in practice. Write (2.4) as 



E[M - t It < M] 
n 1 n 



e nXM - 1 + nXM 
(nX) (l-e'" XM ) 
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and then expand in Taylor’s series to obtain 



P [M — t It < M] 
n 1 n 



(nXM) : 
2 ! 



+ o(nXM) 



nX (nXMfo (nXM) ) 



M 

2 



(2.5) 



This states that in the limit of vanishingly small accident or disease 
rate and average of one-half the mission time will be lost, provided 
that the probability is very small of an accident or breakdown of the 
sort envisioned (e.g. a heart attack, stroke, severe injury, in the 
medical case) . 

Model 2 . Chance Demand, Differing Demand Rates. 

We can realistically generalize Model 1 to the situation for 
which each individual has a different characteristic rate or probability 
of requiring the vital service: A^d t is essentially the demand 

probability for individual j (j = l,2,...,n). Then the overall 
demand rate is, assuming independence, equal to 

A(n) = A + A + ... A (2.6) 

1 Z n 

and the distribution of is still exponential, with A(n) replacing 

nA in (2.2). All this is well-known; see Feller [1]. 

Practically speaking, one might consider classes of individuals 



who are more or less susceptible to disease or accident, and who can 
be characterized by the same failure or catastrophe rates within the 
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class. For example, the older officers (ship captain, executive 
officer, etc.) are probably less prone to incur an appendicitis outbreak 
than would a younger man; on the other hand, an older person might be 
more prone to heart attack. This sort of consideration would compel 
us to put individuals into classes, with characteristic rates 
(A^, u = 1,2, ...,c}; then 

c 

A(n) = l n(u) A ' (2.7) 

u=l u 

where n(u) is the number of individuals in the u — class. 

BuMed data might well be available to provide estimates for 
the above rates. For equipments, 3M data might well be interpreted 
for the same purpose. 

Model 3 . Simple Chance Demand, Demand Rates Randomly Selected from a 
Population. 

An interesting and realistic generalization of Model 2 involves 
the assumption that each individual's "failure" or demand rate is 
itself randomly and independently drawn from some fixed population. 

An immediate generalization, in the spirit of Model 2, is that 
certain classes of individuals 1 demand rates come from different 
populations . 

To explore this sort of assumption we first write down the 
survival probability of the group (crew), given their failure rates. 

This turns out to be, by independence, equal to 
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p {x n > x | \ lf X 2 X n > = e ■* 

Now by assumption each A^ is independently drawn from a population 
with distribution H(y) , and density h(y). Consequently, removal 
of the condition on the A's amounts to integrating: 




P{x > x} 
n 



n 

n 

j=i 



oo 

e -yx h(y)dy = 

0 



[h(x)] n 



where h(x) represents the Laplace transform of the density h. 



Example . Let h be a Gamma density: 



h(y) 



-ay , \ 3~1 

e _ (ay) a 

r (3) 



( 2 . 8 ) 



Then 



■ & 



(2.9) 



and the chance that a mission proceeds with no demand is 



P{t > M} = 
n 




( 2 . 10 ) 



It will also be interesting to derive the expected lost mission time 



under this model. It is 
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M 



E[(M-T n ) + | X 1 ,X 2 



y • • • y 




o 




^ [1 - (l+X (n)M)e X(n)M 




( 2 . 11 ) 



Next the conditions on the X’s must be removed. This can be done 
easily in terms of our gamma density illustration. 

Example , h is Gamma. 

Then A(n) is Gamma with parameters a and Bn. Consequently, 

we need 



00 




0 



oo 



a e 




0 



a 



( 2 . 12 ) 



Bn - 1 



and 
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; E(n) 



-A (n)M 



1 -yM -ay (ay) 

— e e 3 . a dy 

y T(Bn) 3 






Bn 7 x e -y(a+M) 



17 ( 



a+M( a 



3n 



— f — (2 13) 

gn-l^qW U ‘ XJ; 



Assembling the expression for the expectation of (2.11) we find that 



E [ (M-t ) + ] = M 
n 



a 

Bn-1 



+ 



1 

Bn-1 



(q+M) 




(2.14) 



A generalization can be carried out for the case in which several sub- 
groups of crew members are described by their specific gamma distribu- 



tions, but this step will be postponed. 
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3 . Costs and Decisions . 

Armed with various models that describe demand for service by 
a medical doctor or other specialist we can formulate decision analyses. 
Our demand models provide inputs to these analyses, as do certain costs. 
Decision Model 1 . 

Suppose ship missions are of approximately constant duration 
M. Let D be the (dollar) cost per unit time of maintaining a 
medical doctor aboard ship. Then MD is the dollar cost per mission 
of keeping the doctor aboard ship while the ship is engaged in an 
active mission. 

Let R be the cost of the evacuation or rescue operation 
necessary when an emergency arises and no doctor is present. Think 
of R as being an average cost; clearly this cost will vary with the 
location of the ship, and hence the individual, that is the recipient 
of the rescue attempt. 

Let p(n;M) denote the expected number of emergencies that 
arise when a crew of n individuals embark on a mission of duration 
M. Our models of Section 2 provide various bases for evaluating 
p(n;M) . 

Model 1 implies that the number of demands during a mission is 

-XM 

binomially distributed with mean p(n;M) = n(l-e ). Suppose that 
each emergency requires a separate rescue or evacuation operation. 

Then the expected cost of rescues or evacuations in p(n;M)R per mission 
if a doctor is not aboard. Presume that if a doctor is aboard all of 
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these can be avoided, but at cost MD . The optimal decision rule is 
then 

-AM 

Carry a doctor if MD < n(l-e )R 

(3.1) 

-AM 

Do not carry a doctor if MD > n(l-e )R 

If A is quite small, as should often be true, this becomes a good 
approximation! 



Carry a doctor if D < n A R 



Do not carry a doctor if D > n A R. 



(3.2) 



Of course if there is equality (D = nAR, for example) then other 
considerations must settle the matter. 



Decision Model 2 . 

This model simply recognizes the differences between demand 

(injury, accident, or sickness) rates between individuals, as in Demand 

Model 2. For that model the expected number of demands is 

n -A .M 

p(n;M) = l (1-e 2 ) 

j=l 

Hence our decision rule becomes 









n 


-A .M 


Carry a doctor 


if 


HD < 


R l 


(1-e J ). 








j=i 










n 


-X .M 


Do not carry a 


doctor if 


MD > 


R l 


(1-e 2 ). 



j=i 



(3.3) 
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Some additional comments may be made on these models. 

(A) The models tacitly assume that R, the cost of a rescue 
operation, is the same regardless of mission. In fact, one can assign 
a cost that depends upon the mission and then decide on the basis of 
our various decision models whether a doctor can be justified. 

(B) The same comment as in (A) above holds for the rates or 
A -values likely to prevail on different missions. 

(C) The above decision rules, derived for ships, can apply 
also to groups of ships. The doctor can be located on one ship of 
the group, and emergencies will then be transferred to that ship when 
they occur. 

(D) In the above discussion the A-values are taken to be 
known. To make the decision we must obtain estimates, and then treat 
these estimates as equal to the parameter values actually prevailing. 

A more sophisticated approach explicitly recognizes that estimates 
are uncertain; one standard way of handling that situation is by 
means of Bayesian decision theory. We shall apply these ideas in a 
later report. 

The above decision models assume that emergencies generate 
rescue costs, but do not shorten missions. In other situations, perhaps 
having to do with the failure of a major weapon system, this might not 
be the case. It may well be that if a major system goes out on, say, 
a submarine, the latter must return to port prematurely. We set up a 
simple and tentative model for this situation, anticipating that refine- 
ments in the model may suggest themselves. 
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Decision Model 3 . 

Suppose the initial cost for a copy of the ship in question is 
S (dollars) , and that the anticipated life is equivalent to L missions 

S 

of length M. It is reasonable to assess a penalty of — dollars 
per unit time that the ship is not carrying out its assigned task 
during a mission, owing to lack of specialized repair personnel or 
spare parts. 

The expected lost time per mission of length M is obtained 
for Demand Model 1 by multiplying (2.4) by the probability of at least 
one demand during M, namely 1 - e Thus 



Thus the expected cost of lost service if a specialist, or requisite 
spares, are not carried over the life of the ship is 



The optimal decision rule is then derived from the principle that one 
should carry the specialist if his total cost over the life of the 
ship, MLD, is less than the expected cost of curtailed missions: 



-nAM 



E[max(M-T ,0)] = M - 
n 



(3.4) 



Expected cost = 



(lm) L * E[max(M-T n ,0) 




(3.5) 



Carry specialist if 



MLD < S 






nAM 

-nAM 



Do not carry specialist if MLD > S<1 - 



(3.6) 
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A very similar formula can be written down if Demand Model 2 is invoked. 
Note that n now refers to the number of failure— prone equipments 
to be serviced by the specialist. 

One qualitative fact that emerges from (3.6) is that for fixed 
total mission time ML = T one can reduce the need for a specialist 
by shortening mission time, M, (and correspondingly increasing L) . 

By indefinitely shortening M the right-hand side of (3.6) can be 
brought very close to zero, which guarantees that our decision rule 
will recommend that the specialist be left ashore. Of course, indefi- 
nite shortening of the mission time is impractical, but the tendency 
is of interest and can be quantitatively assessed by use of formulas 
like (3.6). 
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